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Abstract. Using the method of separation of variables and a new approach to 
calculations of the prolate spheroidal wave functions, we study the optical properties 
of very elongated (cigar- like) spheroidal particles. 

A comparison of extinction efficiency factors of prolate spheroids and infinitely 
long circular cylinders is made. For the normal and oblique incidence of radiation, 
the efficiency factors for spheroids converge to some limiting values with an increasing 
aspect ratio a/5 provided particles of the same thickness are considered. These values 
are close to, but do not coincide with the factors for infinite cylinders. The relative 
difference between factors for infinite cylinders and elongated spheroids (a/ 6 ^ 5) 
usually does not exceed 20 % if the following approximate relation between the angle 
of incidence a (in degrees) and the particle refractive index m = n + ki takes the place: 
a k, 50|m -l\ + 5 where 1.2 ^ n ^ 2.0 and k ^ 0.1. 

We show that the quasistatic approximation can be well used for very elongated 
optically soft spheroids of large sizes. 
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1. Introduction 

Rapid calculations of light scattering by non-spherical particles are very important in 
many scientific and engineering applications (see discussion in [Q). The simplest 
model of non-spherical particles — an infinitely long circular cylinder is not physically 
reasonable. However, it looks attractive to find cases when this model could be useful 
because the calculations for infinite cylinders are very simple. Therefore, we compare 
the light scattering by elongated spheroids and infinite cylinders. Previously, such a 
comparison was made by Martin and Voshchinnikov 0] for normal incidence of 
radiation (perpendicular to the rotation axis of a particle, a = 90-*). 

Our consideration is based on the solution to the light scattering problem for 
spheroidal particles by the Separation of Variables Method (SVM) (see [0, @] for details). 
A new type of expansions of the prolate wave functions (Jaffe expansions 0) opens a 
possibility to calculate the optical properties of very elongated (cigar-like) particles. 

In this paper, we study the optical properties of prolate homogeneous spheroids 
with large aspect ratios and compare them with those of infinite circular cylinders and 
of spheroids in quasistatic approximation. 

2. Method 

2.1. Prolate spheroid 

Prolate spheroid is obtained by an ellipse rotation around its major axis. The shape 
of particle is characterized by the aspect ratio a/b where a and b are the major 
and minor semiaxes. SVM's solutions to the electromagnetic problem for spheroids 
published by Asano and Yamamoto § and Farafonov p (see also ^) are fundamentally 
distinguished. In the former, the authors used the Debye potentials to present the 
electromagnetic fields. This approach is similar to the Mie solution for spheres. 
Farafonov chose special combinations of the Debye and Hertz potentials, i.e. the 
potentials used in the solutions for spheres and infinitely long cylinders, respectively. 
All the electromagnetic fields were divided into two parts: the axisymmetric part 
not depending on the azimuthal angle (p, and the non-axisymmetric part when the 
integration over cp gives zero. 

The computational efficiency of the Asano and Yamamoto's and Farafonov's 
solutions can be compared in the following way (see |^ for details). Let N is the number 
of terms in sums which give the efficiency factors for extinction and scattering Qext, Qsca. 
with some accuracy in the Farafonov's solution. Then, to obtain the results with the 
same accuracy, the solution of Asano and Yamamoto requires ~ 2A^ and ~ 5A^ terms 
if the aspect ratio a/b ^ 2 and a/b ~ 10, correspondingly. Because the computational 
time is proportional to in this case and t (x N"^ for the Farafonov's solution, the 
advantage of the latter is evident, especially for large aspect ratios. 

It should also be noted that the convergence of Farafonov's solution for spheroids 
follows that of the Mie solution for spheres (see Table 2 in Voshchinnikov, 0). 
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It must be emphasized that the previous calculations for spheroidal particles with 
the Asano and Yamamoto's solution were restricted by the aspect ratios a/b < 5 (size 
parameter 27ra/A < 30) ^ or a/b = 10 {2na/X ^10) fill. 



The particle size can also be specified by the parameters xy = 27rry/A {ry is the 

radius of a sphere of volume equal to that of spheroid, A the wavelength of incident 

radiation) or c = 27i/X ■ (i/2 {d is the focal distance of a spheroid). The expressions 

relating different parameters have the form 

2vra ^ , ^ 

^ = ^0 (1) 

2TTrv 2Tia fa\-'^/^ , , 

where Vy = ab^. The parameter depends only on the aspect ratio a/b 

(3) 
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All cross-sections C (extinction, scattering etc) are connected with the 
corresponding efficiency factors Q via the relation 

C = GQ (4) 

where G is the "viewing" geometrical cross-section of a particle (the area of the particle 
shadow). The factors Q and geometrical cross-sections G depend on the angle of 
incidence a (the angle between the direction of radiation incidence and the particle 
rotation axis, 0*^ < a < 90°). The geometrical cross-section of prolate spheroid is 

G{a) = nb (a^ sin^ a + b'^ cos^ ^ . (5) 

In the case of oblique radiation incidence (a ^ 0°), the factors depend on the state 
of polarization of incident radiation: Q™'"""^. The superscript TM (TE) is related to 
the case when the electric vector E of the incident radiation is parallel (perpendicular) 
to the plane defined by the rotation axis of a spheroid and the wave vector. Expressions 
for the factors can be found in 0. 

In order to compare the optical properties of the particles of different shapes it is 
convenient to consider the ratios of the cross-sections for spheroids to the geometrical 
cross-sections of the equal volume spheres, C/nry. They can be found as 

C [(a/6)2sin^a + cos^a]^/^ 



7rr2 (a/6)2/3 
2.2. Infinite cylinder 



Q. (6) 



The extinction, scattering and absorption cross-sections of an infinite cylinder are, of 
course, infinite. In order to keep some physical sense the infinite cylinder is replaced by 
a very-very long finite cylinder and the cross-sections are calculated per unit length of 
an infinite cylinder (see |T|, p. 203). 



Light scattering by an elongated particle 



4 



For the infinite circular cylinder, the efficiencies are usually defined as the 
corresponding cross-sections per unit length, Ccyi, divided by 2acyi, the diameter of the 



cylinder [12 



For a finite cylinder of length 2L and radius Ocyi, Q is the ratio of the cross-section 
Ccyi to the normally projected geometric area of the cylinder, 2acyi2L. 

The extinction efficiencies are considered for two cases of polarization of incident 



radiation: E case (TM mode) and H case (TE mode) [l|], ||T2 



qE^ ^ ^cyi ^ ^cyi ^ Re IbQ + 2^b^ \ (7) 
2cicyi 2acyi2L Xcyi [ 



2cicyi 2acyi2L Xcyi \ ^ 

where Xcyi = 27racyi/A. The coefficients and 6f depend on the complex refractive 
index m = n + ki, size parameter Xcyi and the angle a (0° < a < 90°). The expressions 



for coefficients are given in |12 



Because of different definitions of geometrical cross-sections for spheroids and 
cylinders, the comparison of factors requires the normalization 

^TM,TE _ Qcy\ /qx 

where a 7^ 0°. 

3. Numerical results 

In previous attempts to find the limits of applications of the model of infinite cylinders 
(0' d])' particles of the same thickness were considered, i.e. spheroids and cylinders 
had the equal parameters 2'nb/X = Xcyi = 27racyi/A. Martin notes that the extinction 
factors for spheroids become to resemble those for cylinders if a/b ^ 4 for the normal 
incidence of radiation {a = 90°). But in order to align the peaks in extinction, the x 
scale for cylinders was stretched by a factor 1.13. 

We suggest another way for comparison of the efficiency factors: to put equal the 
volume and the aspect ratio of a spheroid and a very long cylinder, respectively, i.e. 

Kph = Kyi, T = — • (10) 

Ocyl 

From equation ([T0|) , we obtain the relation between the size parameters for infinite 
cylinder and prolate spheroid 

27racyi _ 27r6 /2\^^^ 

X ~~ [3) ■ 

Note that the scaling factor arising in equation (|Il|) [(3/2)^/^ ^ 1.145] is close to that 
found empirically in 0]. 
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3.1. Normal and oblique incidence of radiation 

The extinction efficiency factors for spheroids and infinite cyhnders are plotted in 
figure |l] for the case of normal incidence of radiation (perpendicular to the rotation 
axis, a = 90°). There is a similarity of the behaviour of the efficiency factors: both 
TM and TE modes converge to some limiting values which are close to those of infinite 
cylinders but do not coincide with them. This occurs if the size parameter is defined 




Figure 1. Extinction efficiency factors (TM and TE modes) for prolate spheroids 
{Qext = Cext/G{a)) and infinite cylinders (Qoxt = C'cxt/(2acyisina)) as a function of 
the size parameter 27racyi/A = 27r6/A • (2/3)^/^. 

by the relation ( [TT| ) and the aspect ratio a/6 ^ 5 (a noticeable difference in factors for 
a/6 = 2 it is seen in lower panels of Figs. |l|, |^). 

If we reduce the angle a (figure ^ upper and middle panels), the similarity in 
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27Taeyi/X=27Tb /A*(2/3) 



Figure 2. Extinction efficiency factors (TM mode) for prolate spheroids (Qext = 
Ccxt/G{a)) and infinite cylinders (Qext = C'ext/(2acyisinQ)) as a function of the size 
parameter 27racyi/A = 27r6/A • (2/3)^/^. 
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behaviour of factors remains although the ripple-hke structure changes. The hkeness 
disappears if we approach the grazing case: the first peak for cyhnders occurs at smaller 
size parameters than for spheroids (see also figure^). For refractive index m = 1.5 + O.Oi, 
it takes place if the angle a becomes smaller than ~ 30°. Because the ripples reduce with 
the decrease of the real part of refractive index n or the increase of its imaginary part k, 
in these cases the infinite cylinders can be used instead of elongated spheroids at smaller 
values of a. However, the difference becomes very large for a ^ 15°. If 1.2 ^ n ^ 2.0 
and k ^ 0.1, the linear relation between the angle of incidence a (in degrees) and the 
particle refractive index m = n + ki takes the place: a ^ 50|m — 1| + 5. Then the relative 
discrepancy between factors for infinite cylinders and elongated spheroids (a/6 ^ 5) does 
not exceed 20 % near the first maximum. 

Anyway, we can expect the resemblance of the extinction properties of very 
elongated particles (spheroids and cylinders) in the case of radiation incidence close 
to normal if the size parameters are recalculated according to equation (p!T|). 

3.2. Parallel incidence of radiation 

Figure ^ shows the extinction efficiency factors for very long spheroids in the case of 
the incidence of radiation along the rotation axis of the particles. The behaviour of 
the factors is rather regular, and their values smoothly reduce with the size parameter. 
There are 25 maxima on the interval 27ra/A = — 300 that is totally distinct from the 

500 




50 100 150 200 250 300 



2TTa/A 

Figure 3. Extinction efRciency factors for elongated spheroidal particles in 
dependence of size parameter 27ra/A. 

extinction by spherical particles (cf. figure ^). The positions of the maxima for large 
particles are determined by the path of light inside the spheroids (the phase shift of the 
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Xv=2TTrv/A 

Figure 4. Normalized extinction cross-sections in dependence of the size parameter 
27rrv / A for spheres and prolate spheroids. 



central ray) . Note that the "equivolume" size of particles xy considered on figure ^ is 
moderate: from equation (0) follows that xy ~ 13.92 if 27ra/A = 300 and a/b = 100. 
But the path of rays inside spheroidal particle is in 300/13.92 ^ 21.5 times longer in 
comparison with spherical particle. 

The values of factors on figure |^ are very large. This is the result of normalization 
by geometrical cross-section which is small in this case and is equal to G{0^) = vrfe^ (see 
equation (|^)). Indicate that the limiting value of extinction factors for particles of any 
shape must be equal to 2 ( "extinction paradox" ) but as it is seen from the insertion on 
figure 1^, this condition is far to be satisfied although the tendency for reduction of Qext 
is observed. 

Figure § shows the normalized extinction cross-sections for spheres and prolate 
spheroids at parallel incidence of radiation [a = 0^). As follows from figure § spheres 
of the same volume scatter more radiation than spheroids but the situation changes 
in the case of another orientation of spheroids. Note also on the absence of ripple-like 
structure on curves plotted for spheroids. 

The ripple-like structure disappear on the extinction curves for cylinders as well 
in the case of nearly grazing incidence of radiation (figure |^). From figure ^, it is 
seen that the similarity between elongated spheroids and infinite cylinders is absent. 
For tangentially incident radiation, the positions and strengths of maxima for infinite 
cylinders do not coincide with those for spheroids at a = 0°. 
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27Taeyi/A = 2TTb/A*(2/3) 



Figure 5. Extinction efficiency factors (TM mode) for prolate spheroids (Qext = 
Cext/G(a)) and infinite cylinders (Qoxt = C'ext/(2acyisina)) as a function of the size 
parameter 2nacyi/\ = 27r6/A ■ (2/3)^/^. 



3.3. Quasistatic approximation 

The optical properties of extremely prolate and oblate particles may be approximately 
calculated using the quasistatic approximation. This is a generalization of the Rayleigh 
and Rayleigh-Gans approximations when the electromagnetic field inside a particle is 
replaced by the incident radiation field (as in the Rayleigh-Gans approximation), taking 
into account the polarizability of the particle (as in the Rayleigh approximation). The 
expressions for the efficiency factors and amplitude matrices were obtained in |T3[. The 
range of validity of the quasistatic approximation is discussed in [Q. 

Figures |] and ^ show the extinction efficiency factors for spheroids calculated with 
the exact solution and the quasistatic approximation in the case of normal {a = 90°) 
and parallel {a = 0°) incident radiation. As we consider non-absorbing particles, the 
extinction efficiencies coincide with the scattering ones which are in the quasistatic 
approximation: 

gTM^QO) ^ c%^(^g-l) |-^|2 r(sinV + cos2^cosV)G'2(w)sinMM99 (12) 
Qtt jo Jo 

Q™(90°) = '^^f^\as\' r r G^{u) sm^edOd^ (13) 
yvr Jo Jo 

where ai and are the polarizabilities, 

3 

G{u) = —{sinu — ucosu) (14) 
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Figure 6. The extinction efficiencies for prolate spheroids and infinite cylinders. 
The results for spheroids were calculated using exact solution and quasistatic 
approximation. 



and 



u = c^o\Kcos9-iy+{^)\m^e (15) 



if a = 0" and 



u = c^f^^cos^e + i^^j (sin2^ + 1 -2sinecos(^) (16) 
if a = 90°. 

Figure ^ shows the enlarged part of figure ^ (lower panel) for small size parameters. 
Here, the results for quasistatic approximation are also plotted. It is seen that this 
approximation describes rather well the behaviour of efficiencies for relatively large size 
parameters (the value 27racyi/A = 1 corresponds to 27ra/A = 5.725 for a/b = 5 and 
27ra/A = 11.45 for a/b = 10). 

The extinction efficiencies for optically soft particles are drawn in figure |^. The 
results given here show the light scattering by a spheroidal cavity (mcavity = 1.49 + O.Oi) 
in glass (mmodium = 1.48 + 0.0^). As follows from figure the quasistatic approximation 
allows to calculate the optical properties with the relative error smaller than ~ 25 % in 
a wide range of particle sizes. 

The range of validity of the quasistatic approximation is discussed in |14|. It was 
found that a maximum parameter xy for which the quasistatic approximation and 
exact theory yielded the results coinciding within 1%, can be described by the following 
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Figure 7. The extinction efficiency factors (left y-axis) for prolate spheroids 

('TT'inedium = 1-48 + O.Oi, rTZcavity = 1-49 + O.Oi, u/b = 10, tt = 0^) in dependence 
of the size parameter 27ra/Amedium- The calculations were made using exact solution 
and quasistatic approximation. The relative error of the calculations of factors in 
quasistatic approximation is given by the right y-axis. 



approximate formulae: 

0.021n(a/6) +0.13 

(n- 1, 

for prolate spheroids and a = 0° 



^ — ^rHiso — (17) 



0.10 

^ _ X)0.131n(a/6)+0.29 ^ °) 



for prolate spheroids and a = 90°, 



< ^iH (19) 

~ (n- 1)0.09 ln(a/6)+0.26 ^ ^ 



for oblate spheroids and a = 0°, 



^ 0.06 ln(a/6) +0.12 

^ (^ _ 1)0.23 (20) 



for oblate spheroids and a — 90°. 



4. Conclusions 

We applied a new method of calculations of the spheroidal wave functions to the study 
of the optical properties of very elongated (cigar-like) spheroids. New approach allowed 
to compare light scattering by prolate spheroids and infinitely long circular cylinders 
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and investigate the applicability of the quasistatic approximation for very elongated 
spheroids. 

It is found that the efficiency factors for spheroids and cylinders have quite similar 
behaviour for the normal and oblique incidence of radiation, if the aspect ratio of 
spheroids a/b ^ 5. The resemblance of factors arises provided spheroids and very long 
cylinders of the same volume and aspect ratio are considered. The following approximate 
relation between the angle of incidence a (in degrees) and the particle refractive index 
m — n + ki takes the place: a ^ 50|m — 1| + 5 where 1.2 ^ n ^ 2.0 and k ^ 0.1. In 
this case, the relative discrepancy between factors for infinite cyhnders and elongated 
spheroids {a/b ^ 5) does not exceed 20% near the first maximum. 

It is shown that the quasistatic approximation rather well describes the extinction 
by very elongated optically soft spheroids of large sizes. 
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